Relativistic Electron Vortex Beams: Angular Momentum and Spin-Orbit Interaction 
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Motivated by the recent discovery of electron vortex beams carrying orbital angular momentum 
(AM), we construct exact Bessel-beam solutions of the Dirac equation. They describe relativistic 
and nonparaxial corrections to the scalar electron beams. We describe the spin and orbital AM of 
the electron with Berry-phase corrections and predict the intrinsic spin-orbit coupling in free space. 
This can be observed as a spin-dependent probability distribution of the focused electron vortex 
beams. Moreover, the magnetic moment is calculated, which shows different g-factors for spin and 
orbital AM and also contains the Berry-phase correction. 
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Introduction. — Vortex beams carrying quantized or- 
bital AM (0AM) along their axes of propagation are 
widely studied and exploited in modern optics [1]. Some 
of us previously predicted the existence of free-space vor- 
tex beams for non- relativistic scalar electrons [2]. More 
recently, electron vortex beams have been produced ex- 
perimentally by several groups [3-5], using spiral phase 
plates [3] and nanofabricated diffraction holograms [4,5] 
which allow generation of beams with an 0AM up to 
100?i [5]. This opens a promising avenue in electron mi- 
croscopy [4], can be employed in quantum problems in- 
volving the AM interaction with external fields [2] and 
other particles [6] , and deepens the analogy between light 
and matter waves. 

The above studies dealt with the simplest approxi- 
mation of a scalar electron and paraxial wave equation. 
However, a self-consistent description of the AM prop- 
erties of the electron requires the exact Dirac equation, 
which describes the spin AM (SAM) intrinsically, and 
takes relativistic effects into account. Indeed, in optics, 
the full vector Maxwell equations are necessary to de- 
scribe nonparaxial vortex beams and the spin-orbit in- 
teraction (SOI) phenomena arising upon the focusing or 
scattering of vortex beams [7-11]. Although zero-order 
Dirac wave packets are studied [12] and higher-order 
modes are used in confined guiding potentials [13], the 
0AM vortex solutions for free relativistic electrons have 
not been previously considered. 

In this Letter we construct exact Bessel-beam solutions 
representing the AM eigenstates of a free Dirac electron. 
Using the recent unifying description of the AM and SOI 
of electromagnetic waves [11] and the semiclassical theory 
of Dirac wave packets [14-16], we give a self-consistent 
description of the 0AM and SAM properties of the Dirac 
electron, and predict a number of SOI phenomena which 
can be verified experimentally. 

Bessel beams from the Dirac equation. — We start with 
the Dirac equation [17] in units with c = 1: 



where ip (r, t) is the Dirac bispinor, a and /3 are the 4x4 
Dirac matrices, p = —ihdr is the momentum operator, 
and m is the electron mass. The positive-energy momen- 
tum eigenstates (plane waves) of Eq. (1) are [17] 
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Here cr are the Pauli matrices, k = p/p is the 
momentum-direction vector, and w = (a, f3)'^ , w^w = 1, 
is the 2-component spinor characterizing the electron po- 
larization in the rest frame with E = m. The probability 
density and current are defined as 

P = '^H , j = , (4) 

so that Pp — 1 and jp ~ p/E for the plane wave (2) and 
(3). The basic states of the polarization, w, can be chosen 
as eigenstates of tr^, i.e., the spin z-component in the rest 
frame [17]. Below we use such states w^: w^^"^ = (1,0)'^ 
and w^^/"^ — (0, 1)"^, with eigenvalues s = ±1/2. 

Bessel beams [18] represent a superposition of monoen- 
ergetic {p = const) plane waves forming a fixed polar 
angle 6 — Oq with the z-axis. Fig. la. They are uni- 
formly distributed over the azimuthal angle (j) G (0, 27r) 
with a vortex phase dependence e'^"^, £ = 0,±1,±2,.... 
Using cylindrical coordinates (r, (p, z) in real space and 
{p±,(j),p\\) — {p sin 9, (j),p cos 9) in momentum space, we 
write the Fourier spectrum for the Bessel beam as: 

(P) - -i^W{p) 5{p^ - p^o) e''^ , (5) 

where p±Q — psin^o, P||o = pcos9q. The beam field is 
given by the Fourier integral of ipi (p) , which yields 



ihdt'ip = (a • p + /3 m) -0 , 
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where ^ {z,t) = h ^ {p\\o^ — Et^ and ^ = k±Qr, k±o = polarization amplitudes w = {a, (3)^ are the same for all 
p^o/h. Substituting here Eq. (3) and assuming that the the plane waves, we evaluate the integral (6) resulting in 
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FIG. 1: (Color online) (a) The Bessel-beam spectrum (5) 
forms a cone of plane waves with a fixed polar angle 9 — 9o 
and color-coded azimuthal phase e'^** on the monoenergetic 
sphere p = const, (b) An example of the probability density 
and azimuthal current distributions for the scalar Bessel beam 
with £ — 1 (the dashed line shows the zero of the current). 



where A = (^1 - — j sin ^o- 

Equation (7) describes spinor Bessel beams for Dirac 
electrons, i.e., the electron counterpart of optical vector 
Bessel beams [11,19]. Hereafter, we use modes -01 with 
polarizations w = w^, s = ±1/2. The first term in the 
square brackets represents a scalar-like Bessel beam of 
the order of i: ipi oc (^) 6'^^"^+*^ Fig. lb, whereas 
the terms proportional to -s/A describe the polarization- 
dependent coupling to other modes of the order of {£ + 
2s). Thus, the latter terms describe the intrinsic SOI, 
determined by the strength A, which vanishes both in 
the paraxial {9q — >■ 0) and non-relativistic (p — > 0) limits 
as A « Ogp/m. Using the canonical 0AM and SAM 
operators, 



L = r X p , S = — diag(cr, cr) 



(8) 



one can readily see that in the limit A ^ 0, the solutions 
(7) are eigenmodes of both = —ihdip and S^, with 
eigenvalues £ and s, respectively. However, in the general 
nonparaxial relativistic case A ^ 0, the Bessel beams ipf 
are eigenmodes of the total AM Jz = Lz + ^z- Jzi^t = 
{£ + s) V'l, but not of the 0AM and SAM separately, cf. 
[11,19]. This also demonstrates the SOI. 



The intrinsic SOI of a free Dirac electron manifests it- 
self in the spin- dependent probability density and current 
distributions, cf. [11]. Substituting the spinorial Bessel 
beam tp^, Eq. (7), into Eq. (4), we obtain 
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JKO - ^J|(02 + ^^ae)j£+2s(0'^, (10) 

with ip and z being the corresponding unit vectors. The 
distributions pg and are invariant w.r.t. the trans- 
formation {£,s) — > {—£,—s), but neither w.r.t. (^, s) — > 
{£, ~s) nor {£, s) — ?> (— s), which reflects the l-s symme- 
try typical of SOI. Figure 2 shows the fine spin- dependent 
splitting of the density profiles p| (^) originated from the 
second SOI term in Eq. (9). In the important particular 
case \£\ — 1, the SOI term has A radically different be- 
haviour at ^ = for the states 2s = £ and 2s — —£, and 
the effect is drastically enhanced. Then the SOI term is 
proportional to Jq (^) and J| {£,) for the two states, re- 
sulting in, respectively, a finite (despite the vortex!) and 
a zero intensity in the beam center, Fig. 2b. This effect 
can be observed for reasonable experimental parameters, 
akin to the tight focusing of circularly polarized opti- 
cal beams with \£\ = 1 [8,9,20]. The azimuthal current 
Eq. (10), exhibits similar behaviour depending on 
the {£, s) quantum numbers. Fig. 2. 

Semiclassical theory of the electron AM. — The ex- 
pectation values of the 0AM and SAM for the Bessel 
beams (7) can be determined using canonical operators 
(8). However, to reveal general features of the AM of 
the Dirac electrons, it is useful to consider the Foldy- 
Wouthuysen (FW) momentum representation [17,14,15], 
separating the positive- and negative-energy components 
in the Dirac equation. The FW transformation is a uni- 
tary analogue of the Lorentz boost to the rest frame: 

= Ufw (p) V'', 
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Ulw (°^ ■p + (3m)UFw 



Dirac Hamiltonian (1): 
= 13 E, and yields 
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(ui, 0) for the plane waves (3). 
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FIG. 2: (Color online) Spin-dependent distributions of the 
probability density (9) and azimuthal current (10) for the 
nonparaxial Bessel beams with £ = 3 (a) and I — 1 (b). The 
s = ±1/2 spin states are indicated by "+" and "-" signs, and 
the SOI parameter is A = 0.3, which corresponds to p ~ 2.4m 
and do ~ 7r/4. The drastic difference between pt{0) and Pi{0) 
should be experimentally observable, cf. [9,20]. 



Owing to the momentum dependence, calculations of the 
coordinate-dependent quantities in the FW representa- 
tion brings about additional gauge fields induced by the 
transformation (11) [14,15]. Calculating the operators of 
the observables O — {r, L, S} in the FW representation 
and using projection onto the positive-energy subspace 
p"*" (2x2 upper left sector of the matrix operators), 
we determine the effective operators for electron states, 
which exclude the influence of the negative-energy levels: 

= P+ (uI^OUfw) [14,15,21]. 

In this manner, the operator for the electron position 
reads [14-16]: 



_ • • P X f / 



(12) 



Here A. = ip^ (Up^^dpUpw) is the non-Abelian Berry 
connection (gauge- field), which brings about the effective 
space noncommutativity, T?.-'] = ih^eijkJ^^, with the 
corresponding Berry curvature [14-16] 



1 



TO , _ 

-.+ (1 



E 



k{k ■ a) 



(13) 



In a similar way, we derive the 0AM and SAM operators 
which may be written as (cf. [11]): 
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Here S = p~''(S) = her/ 2 is the non-relativistic electron 
spin operator and 
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is the SOI operator describing the AM part effectively 
transferred from the SAM to the 0AM: C + S = L + S. 

The operators 7?., C, and S describe the corresponding 
observables for the positive-energy electron states. Sim- 
ilar operators were discussed by Pryce in [22]. In the 
relativistic limit m/E — ^ 0, Eqs. (12)-(16) are quite sim- 
ilar to the results obtained in [11] for photons, whereas 
the expression (15) is similar to the expectation value of 
the spin of an electron wave packet derived in [16]. Note 
that the operator S yields the helicity operator /c (k. • S) 
in the relativistic limit and the spin operator S in the 
nonrelativistic limit (1 — m/E) — 0. 

Using the above operators, we calculate the normal- 
ized expectation values of observables for the Bessel 



beams: (O) 



o 
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trum (5) in the FW momentum representation reads 

'ij^l^ cx w^6{p±^ —p±o)e^^'^ and the convolution means 

2D integration in momentum space [11]. As a result we 
arrive at [23] (R) = zz, (P) = pupz. 



(L) =n{£ + A s) z , (S) = n (s - A s) z 



(17) 



Thus, the expectation values (17) are obtained from the 
operators (14)-(16) via the simple mapping: L — > Mz, 
S — >■ Ss z, and A — ^ A z. 

Equations (17) evidence a spin-to-orbit AM conversion 
for nonparaxial relativistic electron beams (cf. the AM 
conversion in optics [7-1 1] ) . This is a manifestation of the 
SOI intimately related to the spin-dependent probability 
density (9) (Fig. 2). Indeed, using S hsz in the Berry 
connection (12) and integrating it along the contour of 
the Bessel-beam spectrum in momentum space (Fig. la), 
we obtain the Berry phase $b = ~ fA - dp = 27rAs. 
The {£, s)-dependent radius of the beam is effectively de- 
scribed by the quantization condition for the cylindrical 
caustic: k±oR'^ = £-\-As = {Lz) /h, which arises from the 
total phase 27r€ -t- $b around the beam [11]. Thus, the 
spin-dependent parts of the 0AM (17) and probability 
density (9) both originate from the Berry phase gained 
after traversing the circular beam spectrum. 

Magnetic moment. — The interaction of localized elec- 
tron states with an external magnetic field is character- 
ized by the magnetic moment originating from the AM. 
The one-electron magnetic moment can be calculated as 



{M)^y {rxj)dV^ I 



pdV 



(18) 



where e < is the electron charge. Using the definition of 
the probability current (4), Eq. (18) can be represented 
as the expectation value (M) = (^/^j M / V') of the 
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operator M = |r x a in the coordinate representation. 
(This operator resembles the 0AM (8) but with the ve- 
locity operator a instead of the momentum p.) Using 
this analogy, one can calculate the effective magnetic- 
moment operator in the FW momentum representation 

for positive-energy electron: Ai = p+ (^Up^^MUFW^ ■ 
Cumbersome but straightforward calculations yield 

M^^iC + 2S). (19) 

Hence, using Eqs. (17), the magnetic moment for the 
Bessel-beam state is 

{M) = ^{£ + 2s-As)z . (20) 

Thus, in complete agreement with our suggestion in [2], 
the magnetic moment is characterized by a g = 1 factor 
for the 0AM and a g = 2 factor for the SAM of electron 
vortex states. At the same time, the SOI Berry-phase 
correction in Eq. (20) disagrees with the previous wave- 
packet calculations [15,16] for the case £ = 0. This is 
explained by the fact that previous approaches neglected 
(L), which does not vanish even at £ = due to non- 
paraxial corrections. The magnetic moment (20) is an 
important observable that manifests itself in the Zeeman 
interaction with an external magnetic field. 

Discussion. — The above operator calculations of the 
expectation values imply a one-electron state, i.e., in fact, 
a wavepacket state of finite extent in the propagation di- 
rection [2]. For z-independent infinite beams, the number 
of electrons diverges, and one can use the linear densities 
of the characteristics per unit z-length [1,7]. We calcu- 
lated the linear densities of the 0AM, SAM, and mag- 
netic moment per one electron per unit z-length: (L^), 
{Sz), and (Mz)- To regularize the corresponding inte- 
grals over the Bessel-beam cross-section, we modulate the 
transverse electron amplitude by a Gaussian of width a, 
and then consider the limit a — >■ oo. Analytic calculations 
for such electron Bessel-Gauss beams [24] give: 

(Lz) ^h{£ + As) , {Sz)=hs, (Mz) ^£ + 2s . (21) 

This coincides with Eqs. (17) and (20) in the paraxial 
limit but differs by terms of order A. Comparative anal- 
ysis of the one-particle expectation values and linear den- 
sities in the nonparaxial beams deserves further consid- 
eration and will be published elsewhere. 

To conclude, we have constructed the exact cylindrical 
solutions of the Dirac equation, which exhibit nontrivial 
spin- and vortex-dependent properties. The intrinsic SOI 
phenomena originate from the Berry-phase terms and can 
be observed experimentally. We have calculated the ob- 
servable 0AM, SAM, and magnetic moment for exact 
relativistic electron vortex states. 
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